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Abstract The tension wood (TW) properties of a 70-year- 
old specimen of Acer sieboldianum Miq. were analyzed by 
using the G-fiber model that was proposed in our previous 
report. The roles of the G-layer on the origins of (1) a high 
tensile growth stress, (2) a large longitudinal Young’s 
modulus, and (3) a high longitudinal drying shrinkage in the 
TW xylem are discussed on the basis of the simulations 
using the G-fiber model. The results suggest that the G- 
layer generates a high tensile stress in the longitudinal di¬ 
rection during xylem maturation; the longitudinal Young’s 
modulus of the green G-layer becomes significantly higher 
than that of the lignified layer; furthermore, the G-layer 
tends to shrink extraordinarily more than that of the 
lignified layer during moisture desorption. 

Key words Cellulose microfibril • Gelatinous fiber • Growth 
stress * Wood cell wall • Reaction wood 


Introduction 

Tension wood (TW) shows abnormal xylem properties 
when compared with the normal wood (NW), e.g., large 
tensile growth stress, high longitudinal Young’s modulus, 
and large longitudinal shrinkage after drying. Some re¬ 
searchers attribute such behaviors to the physical properties 
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of the gelatinous layer (G-layer) through comparing xylem 
properties and anatomical features in the TW xylem. 14 

To verify these ideas, it is required to observe the behav¬ 
iors or physical properties of the G-layer by isolating it 
directly from the lignified layer. However, it is almost im¬ 
possible to obtain the G-layer cylinder without inflicting any 
damage. No matter how a G-layer cylinder is isolated, it is 
still difficult to provide an accurate measurement because 
the G-layer cylinder is too small to be analyzed by a conven¬ 
tional mechanical testing machine. 

The authors consider that simulation using a mathemati¬ 
cal model of the multilayered wood fiber gives one of 
the most effective approaches for estimating the behavior of 
each cell wall constituent as it exists in the cell wall. 5 In our 
previous report, we proposed a structural model of the 
G-fiber consisting of four-layered cylinders (CML + SI + 
S2 + G), and formulated the mechanical behaviors of 
the G-fiber model on the basis of the reinforced-matrix 
hypothesis. 6 

The formula derived in the previous report contains sev¬ 
eral parameters. We need to optimize those parameters to 
obtain a reasonable result when we simulate the observed 
phenomena on the basis of the G-fiber model. Conversely, 
it can be considered that the optimized values of the param¬ 
eters should reflect certain internal properties and fine com¬ 
posite structures of each constituent material in the G-layer. 
In this report, based on the simulation using the G-fiber 
model, we analyze the observed results for the physical 
properties of the TW xylem from an inclined stem of a 70- 
year-old kohauchiwakaede (Acer sieboldianum Miq.), and 
we explain the role of the G-layer on the origin of distinc¬ 
tive xylem properties in the TW. 


Experiment 

Material and method 

A 70-year-old kohauchiwakaede (Acer sieboldianum Miq.), 
grown on the steep slope of a private mountain in Kiyomi- 
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Table 1 . Observed data for a 70-year-old kohauchiwakaede (Acer sieboldiunum Miq.) 


Measuring positions 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Area composition of domain of each tissue (%) measured at low magnification 







Vessel [V] 

7.2 

6.2 

7.3 

7.7 

8.8 

8.0 

9.0 

7.9 

5.1 

9.9 

Ray [R] 

13.0 

17.8 

12.8 

12.6 

12.3 

17.6 

14.7 

12.3 

16.6 

10.8 

Fiber [F] 

79.8 

76.0 

79.9 

79.7 

78.9 

74.4 

76.3 

79.8 

78.3 

79.3 

Area ratios of the wall in the wood fiber domain (%) measured at high magnification 






G-layer [g] 

8.77 

16.4 

7.88 

10.20 

0.507 

1.61 

0 

0 

0 

0 

Lignified layer [s] 

70.2 

60.5 

63.2 

62.5 

69.1 

63.8 

60.6 

62.0 

67.9 

65.9 

Frequency of the fiber cell in the wood fiber domain (/ 0.01 mm 2 ) 








Total wood fiber 

99 

107 

98 

102 

85 

95 

120 

84 

86 

82 

G-fiber 

51 

73 

34 

44 

3 

11 

0 

0 

0 

0 

Relative frequency of G-fiber [</>] 

0.51 

0.68 

0.35 

0.43 

0.04 

0.12 

0 

0 

0 

0 

Physical properties 











Released strain [e*] (%) 

-0.218 

-0.204 

-0.094 

-0.190 

-0.024 

-0.038 

-0.003 

-0.039 

0.004 

-0.020 

Oven-dried shrinkage [a*] (%) 

1.16 

1.10 

0.48 

0.93 

0.24 

0.52 

0.51 

0.50 

0.22 

0.018 

Green Young’s modulus [E^\ (GPa) 

10.78 

9.21 

8.42 

11.07 

8.92 

6.57 

7.15 

6.37 

5.59 

7.35 

Microfibril angle in the S2 layer (degrees) 

Tension wood region 



Normal wood region 



Number of the measured fiber 

32 





32 





Average [6] (degrees) 

23.1 





27.2 





Standard deviation (degrees) 

(2.18) 





(3.38) 






cho, Gifu prefecture, Japan, 14 cm in diameter at breast 
height (DBH), having a leaning stem, was used for the 
experiment. At breast height, ten measuring points were set 
peripherally on the xylem surface of the leaning stem. At 
each point, the released strain of the longitudinal growth 
stress on the xylem surface (e*) was measured using the 
ordinary strain-gauge method in early April 1988. There¬ 
after, rectangular portions, 70 X 10 X 5 mm and 50 X 10 X 
5mm in the longitudinal (L), tangential (T), and radial (R) 
directions, respectively, were sampled away from the upper 
or lower positions at each measuring point of the released 
strain. Then, respective portions were used to obtain the 
tensile Young’s modulus under the green condition ( E *) 
and the longitudinal shrinkage (a*) from green to the oven- 
dried condition. 3 

After that, a transverse section, 10/mi in thickness, was 
cut from each measuring point of the released strain by a 
sliding microtome, and the section was stained by safranin 
and ferric hematoxylin. The section was then mounted on a 
slide glass with the jelly-like compound of gelatin, glycerin, 
and water. By using a light microscope connected to an 
image processor, microscopic images at large and small 
magnification were photographed within the outermost an¬ 
nual ring of the mounted section. From the images at small 
magnification, the area composition domain of each tissue, 
e.g., vessel element (V), ray tissue (R), and wood fiber (F), 
was computed. From the images at large magnification, the 
area ratios of the lignified layer (s), the G-layer (g), and the 
cell lumen in the domain of the wood fiber were deter¬ 
mined. Frequency of the G-fiber per unit area (N g ) in the 
domain of the wood fiber and that of the normal wood fiber 
(N-fiber) (N n ) were also counted. 

Flat-sliced samples, 5 X 5 X 0.015 mm in L, T, and R 
respective directions, were cut from the outermost annual 
rings of both the NW xylem and the highly developed TW 
xylem. Sampled specimens from the TW xylem were 
quickly dried with ethanol, and were ultra-sonicated in wa¬ 
ter to remove the G-layers from the lignified layer. 7 There¬ 


after, the microfibril angles in the middle layer of the sec¬ 
ondary wall (MFA) were measured by the iodine-staining 
method. 8 

Observed results 

Results are summarized in Table 1, which was already 
reported in our previous report. 3 From Table 1, it can be 
clearly understood that the TW xylem shows quite distinc¬ 
tive properties when compared with the NW xylem. It is 
considered that either of the G-layer formation or the rela¬ 
tively small MFA in the S2 layer of the G-fiber would cause 
the distinctive xylem properties in the TW. However, it is 
still unknown which possibility is more concerned with the 
origin of the TW properties, or whether there is some other 
factor that causes the TW properties. In the present article, 
we attempt to answer this question through simulating the 
mechanical behaviors of the G-fiber on the basis of the 
formula derived in our previous report. 6 


Simulation 

G-fiber model 

A schematic model of the typical G-fiber, consisting of the 
compound middle lamella (CML), the SI and S2 layers, and 
the G-layer, are shown in Fig. I. 6 

Parameters in the basic formula 

In this report, we focus on three biomechanical processes in 
the TW xylem: (1) cell wall maturation, (2) elastic defor¬ 
mation due to action of an axial traction under a steady 
moisture condition, and (3) moisture adsorption. The 
G-fiber tends to shrink or expand in its longitudinal or 
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Fig. la-c. Multilayered structure of the gelatinous fiber, a Microscopic 
structure, b a mechanical model, c crosscut surface of the mechanical 
model. Each consists of compound middle lamella ( CML ), the outer¬ 
most layer of the secondary wall (SI), its middle layer (S2), and gelati¬ 
nous layer (G) 

transverse directions when the biomechanical processes 
occur. We denoted the strains of the dimensional changes 
of the single G-fiber in the longitudinal and the diametral 
directions as s L and e T , respectively, which were simulated 
by the formula derived in our previous report. Correctly 
speaking, it is not a model for the behaviors of an isolated 
fiber, because the constitutive equations used in the formu¬ 
lation consider the conditions of shear restraint imposed by 
neighboring fibers. A basic formula to calculate e L and s T 
contains several parameters, which can be categorized into 
groups as follows. 

Anatomical factors 

Anatomical factors are defined as follows: 

r 0 , r 1? r 2 , r 3 : Outer radii in the CML, the SI and S2 layers, and 
the G-layer, respectively. 
r 4 : Innermost radius in the G-fiber. 

p 0 , p 1? p 2 , p 3 : Ratios of the outer radii to the inner radii in the 
CML, the SI and S2 layers, and the G-layer, respectively. 

Po = '•(Ai, Pi = r x lr 2 , p 2 = r 2 lr 3 , p } = r 3 lr 4 . 

h : Thickness of the CML (= r 0 — rf). 

0: Microfibril angle in the S2 layer (MFA). 

Mechanical factors 

Mechanical factors are defined as follows: 

E t , E 2 , E 3 : Young’s moduli of the framework bundles of the 
polysaccharide oriented in the direction parallel to the 
molecular chain of the cellulose in the SI and S2 layers, 
and the G-layer, respectively. 

S l9 S 2 , S 3 : Double shear moduli of the isotropic skeletons of 
the matrix substances in the SI and S2 layers, and the 
G-layer, respectively. 

S Q : Double shear modulus of the CML. 

Internal expansive terms 

Internal expansive terms are defined as follows: 


e[, e 2 , e 3 : Internal strains in the polysaccharide framework 
bundles in the directions parallel to the cellulose molecular 
chains in the SI and S2 layers, and the G-layer, respectively. 
sf, ef, ef: Internal strains in the matrix skeletons in the SI 
and S2 layers, and the G-layer, respectively. Internal strains 
are caused by the changes of the physical state in the cell 
wall. 

Basic equations to calculate the dimensional changes of 
the single G-fiber 

The basic equations that give e L and e T were derived in our 
previous report. 6 

A =/ll(p)£i m +/l2(p)f2 m +/l3(p)£3 m +/l4(p)fl f +/1 5 (P)£2 

+ fi 6 (py 3 +fn{v)P L 

( 1 ) 

f T (= f,| r =rl) =/2l(p)er +/22(p)£ 2 m +/23(p)«3 m + 

/2 5 (P)£2 +/26(P)£3+/2 7 (P)A. 

where a dot over a character represents the derivative by 
an elapsed time t. Coefficients/ n ,/ 12 ,... ,/ 27 are functions of 
p, and p is a parameter vector whose components are p 0 , p 1? 
p 2 , p 3 , 0 , E t , E 2 , E 3 , S 0 , S l9 S 2 , and S 3 . Some parameters 
depend on t during the cell wall maturation, or the moisture 
adsorption. P L represents an axial traction that acts on 
both ends of the G-fiber. We can calculate the dimensional 
change of the single wood fiber by integrating the dif¬ 
ferential Eqs. 1 along the physical state change of the cell 
wall. 

Time-dependent or moisture-dependent behaviors 
of parameters 

Maturation process of the cell wall 

The amorphous constituents, such as xylan and lignin, are 
irreversibly accumulated among the gaps of the polysaccha¬ 
ride bundle after the completion of the polysaccharide 
framework of the CMF and other oriented polyose. In this 
process, the amorphous constituent hardens into the matrix 
skeleton. Thus, S u S 2 , and S 3 tend to increase monotonously 
from very small values to their final values. Moreover, 
the amount of the substance increases irreversibly inside 
the matrix skeleton whose volume is spatially limited. 
As the inevitable consequence, internal strains e 2 , and 
sf are induced in the SI and S2 layers, and the G-layer, 
respectively. 

It is considered that time-dependent changes in E l9 E 2 , 
and E 3 are somewhat smaller than in S t , S 2 , and S 3 because 
the polysaccharide framework had been already completed 
before the matrix substance began to deposit. However, we 
should not ignore the possibility that aging effects, such as 
an increase in the crystallinity of the CMF, would generate 
internal strains e{, s 2 , and s f 3 in the polysaccharide frame¬ 
work bundle. 3,6,911 In such a case, we need to assume a 
certain value for each of them. 
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Drying process in the cell wall 

Because the completed xylem (i.e., green wood) contains 
much water, drying is required before converting the wood 
into natural resources for building or for producing furni¬ 
ture. In this process, the water molecule is discharged from 
the absorption site in the matrix skeleton, after which the 
matrix skeleton tends to shrink and harden. This means that 
S l9 S 2 , S 3 , ef 9 e 2 , and ef to change their values monoto¬ 
nously in accordance with the moisture desorption. At the 
same time, a certain physicochemical change may occur in 
the bundle of the CMF. However, it is natural to consider 
that changes of E l9 E 2 , E 3 and values of e[ 9 s 2 , and £3 are 
somewhat smaller than those of S l9 S 2 , S 3 , sf, t: 2 \ and ef 
because the crystal domain, which is a main component of 
the polysaccharide framework, almost does not participate 
in the adsorption of the water molecule. 

Determination of the values to be assumed for 
parameters in Eqs. 1 


6 is one of the anatomical factors in p. In the present 
simulation, we used the measured values of the MFA in 
the S2 layer of the N- and G-fibers, which are displayed in 
Table 1. 

Mechanical factors 

The SI and S2 layers, and the G-layer can be regarded as 
the parallel composites of the crystalline bundle of cellulose 
and the matrix skeleton. It follows that the simple mixture 
law can be applied to calculate the values of E l9 E 2 , E 3 , S l9 S 2 , 
and S 3 as follows : 12 

E\ — A 1 X C 1 X E cry , E 2 — A 2 X C 2 X E cry , 

E 3 — A 3 X C 3 X E cry , 

c _ (l — A x C^)E maix _ (l — A 2 C 2 )E matr /o\ 

Y 7“7 > Y ^“7 3 ' ' 

1 + V 1 + V 

c _ (l — A 3 C 3 ^E matY 


Anatomical factors p 0 , p h p 2 , p 3 , and 0 


To determine the values of p 0 , p l9 p 2 , and p 3 , it is required to 
know the ratio of the area of each layer to the total crosscut 
area of a single wood fiber. Then, we interrelate the param¬ 
eters p 0 , p l9 p 2 , and p 3 by using the following formulus (see 
Appendix A). 


P0P1P2 





1 


1 - ^ 


Pi( 1 + h/r^l - s ’ ' ]j(l ~ s) - f-g/N s 

(for N g ^ 0), or p 3 = 1 (for N g = 0), 


( 2 ) 


where s and g represent the area ratios of the lignified layer 
(= CML + SI + S2) and the G-layer in the domain of the 
wood fiber, respectively. / and N g represent the numbers of 
the wood fiber and the G-fiber per unit area in the domain 
of the wood fiber. These are experimentally determined 
values. To determine the values of p 0 , p l9 and p 2 by using 
Eqs. 2, we need to give at least two of them. In the present 
calculation, and with reference to previous studies , 5,12 we 
hypothesized 0.025 as the value of h/r l9 and 1.1 as the value 
of p x . Thereafter, for each measuring point of the released 
strain, we calculated the values of p 2 and p 3 by using Eqs. 2. 
The estimated values of p 2 and p 3 are displayed in Table 2. 


where v is Poisson ratio, which is hypothesized to be 0.5 in 
the same way as in our previous reports. ’ “ C l9 C 2 , and C 3 
are crystallinity indices of the polysaccharide framework in 
the SI and S2 layers, and the G-layer, respectively. A l9 A 2 , 
and A 3 are weight ratios of the polysaccharide framework in 
the respective layers. In this study, the values of A l9 A 2 , and 
A 3 are assumed as given in Tablse 3. 

E matr is the Young’s modulus of the molded matrix 
substance, which clearly depends on the elapsed time dur¬ 
ing the cell wall maturation (or moisture content during 
moisture adsorption). With reference to Cousins’s experi¬ 
ments , 14,15 it is assumed that E matr = 2GPa in the green 
condition, and E matr = 4-6 GPa in the dried condition. 
On the other hand, it is considered that the Young’s modu¬ 
lus of the cellulose crystal along the direction parallel to the 
molecular chain (£ cry ) is not affected by moisture ad¬ 
sorption. With reference to the study of Nishino et al ., 16 
we assume E cry = 134 GPa regardless of the moisture 
content. 

Then we assumed the values and Fdependent patterns of 
E l9 E 2 , E 3 , S l9 S 2 , and S 3 as displayed in Table 3 on the basis 
of the above-mentioned discussions and the subsidiary con¬ 
ditions described hereafter, provided that the noncrystal¬ 
line region in the framework bundle was regarded as the 
matrix substance from the mechanical point of view. The 
value of S 0 was calculated by the method described in our 
previous report . 12 


Table 2. Values of the parameters p 2 and p 3 which are estimated from Eqs. 2 


Measuring positions 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Mean 

P2 

P3 

1.625 

1.527 

1.411 

1.598 

1.462 

1.616 

1.448 

1.645 

1.596 

1.367 

1.474 

1.274 

1.413 

1.000 b 

1.439 

1.000 b 

1.558 

1.000 b 

1.519 

1.000 b 

1.494 a 

1.505 c 


a No significant difference among the positions (goodness of fit) 
h No G-fiber was formed at positions 7, 8, 9, and 10 
c Mean value for the six positions where G-fiber is formed 















226 


Table 3. Values of parameters used for the simulation. 


a Chemical composition in each layer 


Layer 

Polysaccharide 

Matrix 

Crystallinity of 


framework (%) 

substance (%) 

the framework 

CML 

15 

85 

100 (%) 

SI 

26 

74 

C\ 

S2 

52 

48 

c 2 

G 

90 

10 

c 3 

b Time-dependent changes of the mechanical properties of constituents 

Time a 

E ay E m J (in SI) 

E matr b (in S2) 

E m J{ inG) S 0 b 


t = 

0 ~ 

' T‘ 

134 

Increase 
from 0 to 2 

0 

0 

4 

t = 

T, ■ 

- T 2 c 

134 

2 

Increase 
from 0 to 2 

0 

4 

t = 

t 2 - 

- t 2 c 

134 

2 

2 

Increase 
from 0 to 2 

4 

t = 

r, ■ 

r ~p d 

1 4 

134 

Increase 
from 2 to 4 

Increase 
from 2 to 4 

Increase 
from 2 to 4 

4 


c Geometrical properties of the layers 

Time a 

Po 

Pi 

Pi 

P 3 e 

Qi.% 

o 

l 

o 

1.025 

1.1 

1.494 

1.0 or 1.505 

27.2 or 23.1 

t = F ~ t 2 c 

1.025 

1.1 

1.494 

1.0 or 1.505 

27.2 or 23.1 

t=T 2 ~ 7V 

1.025 

1.1 

1.494 

1.0 or 1.505 

27.2 or 23.1 

t = T 3 ~ r 4 d 

1.025 

1.1 

1.494 

1.0 or 1.505 

27.2 or 23.1 

CML, compound middle lamella 


a Assumed with reference to Subsidiary condition 1 
b Unit: GPa 

c Maturation process after the formation of the polysaccharide framework in the SI, the S2, and 
the G-layers 

d Drying process from the fiber saturation point (FSP) (t = T 3 ), to the oven-dried state (t = T 4 ). 
r = (T 4 - t)/(T 4 - T 3 ) 

e p 3 = 1.0 in the N-hber, p 3 = 1.505 in the G-hber 
f Unit: degrees 

g 6 = 27.2 in the N-fiber, 0 = 23.1 in the G-hber 


Internal expansive terms 

Neither of the values nor cdependent patterns can be mea¬ 
sured for the internal expansive terms e [, s 2 , s f 3 , e 2 , and 
e™. However, we can optimize their values and f-dependent 
patterns so as to obtain a reasonable simulation. 


Results 

Young’s modulus of the green G-layer 
Experimental results 

Matured secondary xylem of kohauchiwakaede consists of 
four tissue domains, i.e., the wood fiber, the vessel element, 
the ray parenchyma, and the axial parenchyma. It is consid¬ 
ered that these tissues are arranged in a row in the direction 
parallel to the axis of wood fiber, and the following formula 
can be used for calculating the longitudinal Young’s modu¬ 
lus of the TW xylem ( E *) by the simple law of mixture. 



- (EEl + V-EZ + REl 

F + V + R + P y 




where E F L , E V L , E \, and E P L are Young’s moduli of the wood 
fiber, vessel element, ray parenchyma, and axial paren¬ 
chyma, respectively, under the green condition, and F + V 
+ R + P = 1. Considering E^/E[ << 1, E R L /E F L << 1, and 
E p l /E f l << 1, we obtain E F L = E X JF. In the case of 
kohauchiwakaede, the amount of the axial parenchyma is 
somewhat less than that of the other tissues, and its mor¬ 
phological features are similar to the wood fiber cell apart 
from the fact that the wall thickness of the axial paren¬ 
chyma is smaller than that of the wood fiber. In this study, 
for simplification, we did not distinguish the axial paren¬ 
chyma from the wood fiber when we determined the values 
of F , V, R, and P. 

According to the observations, there was no significant 
difference among the measuring points on the periphery for 
the morphological properties of the G-fiber, e.g., the thick¬ 
ness of the lignified layer, that of the G-layer, and their 
morphological appearances. The same can be said in the 
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case of the N-fiber. Then, by applying the simple mixture 
law to the fiber domain that is regarded as a parallel com¬ 
posite of the G-fiber and the N-fiber, we obtain the follow¬ 
ing formula: 

El = <p • E[ + (1 - <p)E[ = (E[ - El)<p + El, (5) 

where <p = Njf, N g + N n =f 

El and E n L are the axial Young’s modulus of the green 
G-fiber and that of the green N-fiber, respectively, and 0 is 
the relative frequency of the G-fiber in the fiber domain. On 
the other hand, we obtained the relationship between 0 and 
E f l (= E x l /F) as shown in Table 1, which was approximated 
by the following linear regression: 

El = 7.740 + 8.50 (r = 0.857). (6) 

Then, comparing Eqs 5 and 6 directly, we obtain 

El = 16.24 GPa, El = 8.50 GPa, (7) 

provided that we do not use the data obtained from mea¬ 
suring point 5 when deriving Eq. 6. This is because the 
observed value of the longitudinal Young’s modulus at 
measuring point 5 was larger regardless of having very small 
amount of G-fiber formation, and therefore the estimated 
value of Ef becomes abnormally large at position 5 in com¬ 
parison with that at the other positions. It is reasonable to 
assume that some error occurred in measuring the elastic 
modulus of the specimen at position 5. 


Simulation using the wood fiber model 

In this simulation, we assumed the condition of the steady 
moisture state (green condition, i.e., the state at t = T 3 in 
Table 3). Then every component in p must be constant, and 
both dsf and de\ (i = 1, 2, 3) should all be nil. Then, from 
Eqs. 1, we obtain the following formula to calculate the 
longitudinal Young’s modulus of the wood fiber (E L ): 

e l = {i /(nro)}dP L /de L = {l/ (jr r 0 2 )}// 17 (p). (8) 

The values assumed in Table 3 were used for the simulation 
using Eq. 8. Firstly, we optimized the values of C, and C 2 in 
Eqs 3 to simulate the experimentally determined value of 
El (= 8.50 GPa). In this simulation, we assumed that the 
degree of crystallinity in the framework bundle of the ori¬ 
ented polysaccharide was identical in the SI and the S2 for 
convenience, because there is no reason for considering that 
properties of the CMF are different in the SI and the S2 
layers. Thereafter, we applied the optimized values of C, 
and C 2 to the simulation of E L in the green G-fiber, and 
optimized the values of C 3 to obtain the experimentally 
determined value of Ef (= 16.24GPa). Finally, the opti¬ 
mized values of C 1? C 2 , and C 3 became: 

Cj(= C 2 ) = 0.494, C 3 = 0.221. (9) 

From this result, we calculated the longitudinal Young’s 
modulus of the lignified layer in the N-fiber ( E n N ), that of the 
lignified layer in the G-fiber (El), and that of the G-layer 
(E g G ) as follows: 


In the N-fiber: El = 13.13 GPa 

In the G-fiber: E g = 16.28 GPa, E g = 28.27 GPa. 



Growth strain in the G-layer 
Experimental results 

The wood fiber, vessel element, ray parenchyma, and axial 
parenchyma in the differentiating xylem tend to deform 
during their secondary wall maturation. Thus, the growth 
strain is generated in the maturing xylem. The infinitesimal 
increase in the longitudinal growth strain of the xylem at the 
macroscopic level (e *) can be expressed by the following 
formula by the simple mixture law: 


F • El • del + V • El * del + R ■ E* * del 

+ P-E*-del _ (11) 

F ■ El + V El + R El + P El 


where del, de\, de\, and de p L are infinitesimal increases of 
the longitudinal growth strain in the respective tissues. As¬ 
suming E v l /E f l << 1, E r l /E f l << 1, E p l /E f l << 1, and E + V 
+ R + P = 1, we obtain 



Moreover, we obtain the following formula: 


del [= del) 


_ <p-El-del + (l-cp)-El-del 


4>-El + (l-<p)-El 
= •</> + £ l-del 

[El- El \0 + El 



where del and de n L are respective increments in the longitu¬ 
dinal growth strain of the G-fiber and the N-fiber, respec¬ 
tively. By integrating Eq. 12 along the cell wall maturation, 
we can obtain the growth strain of the newly formed xylem 

(£ X l)- 

In order to integrate Eq. 12, we also need to know the 
changes of E[ and E n h during the process of secondary wall 
maturation. It is considered that deposition of the matrix 
constituents has almost no effect on the increases of Ef and 
El because the stiffness of the matrix substance is some¬ 
what smaller than that of the framework bundle. Therefore, 
it is natural to consider that increases of Ef and E\ are 
caused by a certain qualitative change of the CMF, such as 
further crystallization of cellulose. Unfortunately, it is still 
quite difficult to know the time-dependent change of the 
CMF crystallinity in the cell wall. In the present study, for 
convenience, we assumed that the crystallinity in each layer 
is almost unchanged during cell wall maturation. We then 
hypothesized that the El and E\ become constant through 
the cell wall maturation. 

Growth-stress generation is a biomechanical process 
during maturation (lignification) of the secondary wall. 11,17,18 
Thus, we integrate Eq. 12 along the cell wall maturation in 
the G-fiber. As the result, we obtain the following formula: 
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x / f \ = 4>-El-el+ (1-4,)-El-el 
£l1 ^ <t>-El + (1-0) -El 

= (£ g - f £-^-^)-0 + £ L n -^ (B) 

(El- El)-</, + El 

where 

el= jdel, sl= jdel. 

Maturation Maturation 

process process 

Results from Eq. 7 were used as the values of El and E n L 
in this formula. The observed relationship between 0 and s x L 
(= which is shown in Table 1, was approximated by the 
following curvilinear regression: 

= -0.5554 + °' 6003 . f r = 0.956) (14) 

4) + 1.098 v ’ 

By comparing Eqs. 13 and 14 directly, we obtained the 
growth strains of G-fiber (el) and the N-fiber (si) as follows: 

s[ = -0.2693%, el = -0.0087%. (15) 


Simulation based on the G-fiber model 

We integrated the basic formula of Eq. 1 during the G-fiber 
wall maturation under the assumption of dP L = 0. As initial 
conditions, we adopted s L (t) \ t=0 = 0, s T (t) \ t=Q = 0. Results 
from Eq. 9 were used as the values of C 1? C 2 , and C 3 . Values 
of the parameters assumed in Table 3 were also used for 
the calculation. We then optimized the increments and t- 
dependent patterns of e[, s 2 , s\ 9 sf, s 2 , and sf to obtain the 
results in Eq. 15. However, before integrating Eqs. 1, 
we need to know how the maturation of the G-fiber wall 
proceeds. 

Some scientists clarified the lignification process in the 
secondary wall of the softwood tracheid and the hardwood 
normal fiber . 1920 On the other hand, maturation of the G- 
fiber has remained unclear. Recently, based on the tech¬ 
nique of immuno-TEM observation, Kim et al . 21 discovered 
that the activity of the peroxidase is localized in the second¬ 
ary wall after the completion of the G-layer. This suggests 
that lignification proceeds in the secondary wall after the 
formation of the thick G-layer. With reference to those 
investigations, we assumed the following conditions regard¬ 
ing the maturation of the G-fiber. 

Subsidiary condition 1 : Lignification in the SI layer starts at 
t = 0 after the formation of the frameworks of the cellu¬ 
lose and the other oriented polysaccharide in the second¬ 
ary wall and the G-layer, and ends at t = T t . This is the 
first integration interval. Lignification in the S2 layer 
starts at t = 7\, and ends at t = T 2 . This is the second 
integration interval. In the G-layer, deposition of a cer¬ 
tain matrix substance should proceed, however, no 
lignification occurs. In this study, as the third integration 
interval, the deposition of the matrix substance in the G- 
layer starts at t = T 2 and ends at t = T 3 . Then the G-fiber 
maturation is completed at t = T 3 . S l9 S 2 , and S 3 tend to 


increase monotonously and smoothly from very small 
values to their final values in their respective integration 
intervals. 

We integrate Eqs. 1 as follows: 



integration 



integration integration 

interval interval 


We need to impose certain subsidiary conditions on val¬ 
ues and ^-dependent patterns of e[, s [, e 3 , sf, s 2 , and £ 3 m to 
simulate the observed values of s\ and s\. In the case of 
softwood xylem, the observed relationship between the lon¬ 
gitudinal growth strain and the MFA in the latewood trac¬ 
heid can be simulated by supposing [increment in sf] = 1 %, 
[increment in e™] = 0.5%, and [increment in s[] = [incre¬ 
ment in s f 2 ] = —0.15%. 11 With reference to this result, we 
assumed the following subsidiary conditions. 

Subsidiary condition 2: sf and e 2 take positive values. Each 
of them increases monotonously and smoothly from 0 to 
a certain value (= increment) as lignification proceeds in 
each integration interval . 11 It is natural to consider that 
increments in sf and sf depend on the lignin content in 
the respective layers. This is based on the lignin swelling 
hypothesis. However, we assume sf = 0, since no 
lignification occurs in the G-layer. On the other hand, s[ 
and s 2 take negative values. Each of them tends to change 
monotonously and smoothly from 0 to a certain value 
(= increment) with the maturation in each integration 
interval. This postulates the cellulose tension hypothesis, 
which considers that the CMF framework tends to con¬ 
tract in the direction parallel to the cellulose molecular 
chain with the aging of the CMF . 17 22 

Firstly, we simulated the generation of the growth strain 
of the N-fiber (s n L = —0.0087%) by integrating Eq. V 
under the above subsidiary conditions, and optimized the 
increments in s\ and s 2 to obtain the observed value of 
s\ (= —0.0087%). Thereafter, we tried to simulate the gen¬ 
eration of the growth strain of the G-fiber (s g L = —0.2693%) 
and optimized the increment s 3 . In this simulation, we as¬ 
sumed the following subsidiary condition in addition to the 
above two conditions: 

Subsidiary condition 3: According to the observations by 
using a light microscope or an ultraviolet microscope, 
there is no specific difference in the morphological ap¬ 
pearance between the secondary wall of the N-fiber and 
that of the G-fiber. J From this fact, we assumed that t- 
dependent patterns and increments in each of sf, e 2 , s[, 

r 

and s 2 take identical values between the N-fiber and in 
the G-fiber. 

S l9 s[, and sf are all expressed as monotonously increas¬ 
ing (or decreasing) functions of t in the first integration 
interval. S 2 , s f 2 , and e 2 are also monotonously increasing (or 
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Table 4. Combination on the increments in £, m , el, e\, e 2 which optimize the observed value of e n L 


£f (at t = Ef 

1 

0.8 

0.6 

0.4 

0.2 

0 

- 0.2 

-0.4 

el (at t = T 2 y 

0.5 

0.4 

0.3 

0.2 

0.1 

0 

- 0.1 

- 0.2 

£1 (at t = Tl) h (= £2 at t = T 2 ) 

-0.178 

-0.145 

- 0.112 

-0.078 

-0.045 

- 0.012 

0.021 

0.055 

£t (at t = T 2 ) c 

0.170 

0.135 

0.099 

0.063 

0.028 

-0.0078 

-0.043 

-0.079 

Data given as percentages 









a Assumed arbitrarily under Subsidiary conditions 1 and 2 







Calculated from each pair of el 

and el to give the observed value of e\ 






c Calculated from each combination of el, el, 

„f „f 
£ 1> £ 2 








Table 5. 

Estimated value of the increment in £3 

which gives the observed value of e g L 





el (at t -- 

= Ef 

1 

0.8 

0.6 

0.4 

0.2 

0 

- 0.2 

-0.4 

el (at t -- 

= r 2 y 

0.5 

0.4 

0.3 

0.2 

0.1 

0 

- 0.1 

- 0.2 

e\ (at t = 

■- r 1 ) b (=<4atf= t 2 ) 

-0.178 

-0.145 

- 0.112 

-0.078 

-0.045 

- 0.012 

0.021 

0.055 

4 (at t = 

- r 3 y 

-0.721 

-0.737 

-0.754 

-0.770 

-0.786 

-0.802 

-0.818 

-0.834 

£® (at t = 

- Ef 

0.333 

0.298 

0.263 

0.228 

0.193 

0.158 

0.123 

0.088 


Data given as percentages 

a Assumed arbitrarily under Subsidiary conditions 1, 2, and 3 
b Calculated for each pair of el and so as to obtain the observed value of el 
c Calculated for each combination of el, el. e[, e 2 so as to obtain the observed value of s[ 
d Calculated for each combination of el, el, ej, £ 2 , and £3 


decreasing) functions of t in the second integration interval. 
The same can be said for S 3 , e f 3 , and £ 3 m in the third integra¬ 
tion interval. Each of these monotonously increasing 
(or decreasing) functions can be transformed into functions 
that do not contain explicitly 7\, T 2 , and T 3 by transforming 
the integral variable t into y (= t/Ty, 0 < t < Ti), or £ ( = 
it - T t )l(T 2 - T\); T l <t< T 2 ), or k(= (t - T 2 )I(T 3 - T 2 ); 
T 2 <t < r 3 ). Moreover, we know these variable transforma¬ 
tions alter corresponding integration intervals in Eq. 1' into 
an identical one that is from 0 to 1. Thus, the concrete value 
of Eq. V does not depend on 7\, T 2 , and T 3 . Furthermore, 
we should note that integration of Eq. V is not affected by 
the functional shapes of Cdependent variables if each vari¬ 
able changes its value very smoothly in each integration 
interval. This is quite reasonable because we consider that 
the f-dependent changes of these variables gradually 
proceed with the maturation of the matrix skeleton in the 
respective layers (see Appendix B). 

Thus, we can optimize the value of the increment in e 3 as 
displayed in Table 5 which became quite larger than those 
in e[ and e\ as shown in Table 4. 


Drying shrinkage of the G-layer 
Experimental results 

We can describe the shrinking process of wood as a function 
of the moisture content r that is normalized by the moisture 
content at the fiber saturation point (FSP). We denote 
the longitudinal shrinking process of the wood as a x L (r). 
According to the definition, the longitudinal shrinkage 
a\(r) must satisfy the following boundary condition, 
a*(r)| T=1 = 0. a x L (r)\ t = 0 (= a x L ) means the oven-dried 
shrinkage of the wood. An infinitesimal increase of the 
moisture content (dr) causes an infinitesimal change in 
the shrinkage of the wood (da X L ), which is described in the 
following formula: 


F • El • dal + V • £ l v -dal + R- £ L R • dal 
+ P • E[ • dal 

F • El + V • E l v + R • EZ + P ■ E[ 



where dal , da y L , da\ , and dal stand for infinitesimal 
changes of the longitudinal shrinkage in the respective tis¬ 
sues. El> e\, e\ and El are the respective Young’s 
moduli at the moisture content t. Assuming eV El << 1, 
Ell El « 1, El/El « 1, and F+V+R + P = l,v/e 
obtain 


dal — dal. 

We apply the simple mixture law to the fiber domain con¬ 
sisting of the N-fiber and G-fiber in parallel, then, we obtain 
the following formula: 


dal 


( =da l ) 


<p- El ■ da[ + (l — <fi) ■ El ■ dal 


El■da[ 


El 


4>-El + (l-<p)-El 
dal ) ■ (j> + El ■ dal 


[El-El)-<p + El 



where dal an d dal are infinitesimal changes in the shrink¬ 
age of the G-fiber and that of the N-fiber, respectively, and 
Ei and El are the axial Young’s moduli of the G-fiber and 
N-fiber, respectively. 

We obtain an oven-dried shrinkage of the wood fiber 
domain a F L (= a\(r) \ T=0 ) by integrating Eq. 17 from an 
arbitrary r to FSP (r = 1) and extrapolating r —» 0, provided 
that we need to know the r-dependent patterns of Ex. and 
El in advance. Then, we tentatively expressed El and E l 
as follows: 


El = El ■ |(r), El = El ■ £(r) (18) 

where £(t) and £(r) are monotonously decreasing functions 
for r, and they satisfy §(r) | T=1 = 1, and £(r) | T=1 = 1. El and 
El are constants, which stand for the axial Young’s moduli 
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Table 6. Combinations of the increments of e“, e”, £*, e 2 which optimize the observed value of a\ 

El" (=e 2 m ) at t = r 4 ab 15 12 9 6 3 0 -3 -6 

e\ (= 4 ) at t = r 4 bc -3.345 - 2.605 -1.865 -1.126 - 0.386 0.353 1.094 1.833 

a j d,e 9.017 7.361 5.643 3.861 2.010 0.086 -1.915 -3.997 

Data given as percentages 
a Assumed arbitrarily 

b Positive sign means that each cell wall constituent swells during water sorption 
c Calculated from each pair of e“ and ef so as to optimize the observed value of a n L 
d Calculated from each combination of ef, e{, e\ 

e Positive sign means that the wood fiber shrinks in diameter from the FSP to the oven-dried state 


Table 7. Estimated value of the increments in £3 which gives the observed value of a l 


e ; 11 (=£ 2 m = £ 3 m ) at / = 77’ b 

15 

12 

9 

6 

3 

0 

-3 

-6 

£ f ! ( = £ 2 ) at t= r 4 b ’ c 

-3.345 

-2.605 

-1.865 

-1.126 

-0.386 

0.353 

1.094 

1.833 

4 (at t = r 4 ) bd 

5.078 

4.782 

4.486 

4.190 

3.897 

3.603 

3.304 

3.010 


13.66 

11.14 

8.461 

5.616 

2.588 

-0.641 

-4.089 

-7.784 


Data given as percentages 
a Assumed arbitrarily 

b Positive sign means that each cell wall constituent swells during water sorption 
c Calculated from each pair of and ef so as to optimize the observed value of a n L 
d Calculated from each combination of e 2 , e[, e\ so as to obtain the observed value of al 

e Calculated from each combination of e™, ej 1 , e[, s f 2 , and £3 

f Positive sign means that the wood fiber shrinks in diameter from the FSP to the oven-dried state 


of the green N-fiber and the green G-fiber, respectively. For 
simplification, we assumed £(t) = £(t) for all r, which means 
the decreasing pattern of the longitudinal Young’s modulus 
in the G-fiber is similar to that in the N-fiber. Then, Eq. 15 
becomes 


1. For convenience, we calculated the swelling deformation 
of the wood fiber model £ L (x) by integrating Eqs. 1 from 
x = 0 to r = 1. The relationship between swelling e L (= 
e L (r) I r= i) and oven-dried shrinkage a L (= a L (x) | T=0 ) are 
related to each other by the following formulas: 


da[(= dal) = 


[El • dal ~ El • dal ) • <f> + El • dal 


(El- E[)-<p + El 


(19) 


Under these assumptions, we substituted the results from 
Eq. 7 to El and E n L in Eq. 19. As the initial conditions, 
a l( r ) I T =i = a K r ) I T =i = were required. Thus, Eq. 19 can 
be integrated over moisture content (from an arbitrary x to 
r = 1). We obtain the oven-dried shrinkage of the wood 
fiber domain a F L (= a F L (x) | r=0 ) as the following formula. 



(El-al- El-al)-4> + El-al 
(El- El)- <p + El 



The observed relationship between 0 anda^ which is shown 
in Table 1, was approximated by the following curvilinear 
regression: 


a[ = -2.429 + 2363 . (r = 0.867) 

<p + 1.098 v ’ 



Then, comparing Eqs. 20 and 21 directly, we obtain the 
oven-dried shrinkage of the N-fiber (al) and the G-fiber 
(a[) as follows: 


al = 0.2771%, a\ = 1.3026%. 



Simulation based on the G-fiber model 


Free dimensional change of the single wood fiber due to 
moisture adsorption was simulated on the basis of the con¬ 
ditions assumed in Table 3. Thus, dP L should be null in Eqs. 



The integral interval for calculating e L (= s L (x) | r=1 ) is from 
the oven-dried state (t = T 4 ; x = 0) to the FSP (t = T 3 \ 
x — 1). It is regarded that the increasing moisture content x 
is equivalent to the reciprocal elapsed time t. The results 
from Eq. 9 were used as the values of C 1? C 2 , and C 3 in 
this simulation. Then, we optimized the increments in 
e™, £ 3 m , e[, e f 2 , and e 3 to obtain the observed values of a\ and 

a]l 


Swelling of the softwood tracheid cell wall is mainly 
caused by swelling of the matrix substance, e.g., hemicellu- 
lose, lignin, and noncrystalline cellulose. 12,24-27 Therefore, 
it is quite reasonable to postulate that e™, £™, and £ 3 m take 
positive values with the increase of the moisture content, 
and increase monotonously from 0 to the final values, that is 
to say, increments. 

Firstly, we simulated the swelling of the N-fiber 
(e n L = 0.2779%). In doing so, we optimized the values of 
increments in e™, £ 2 m , £[, and e 2 to give the observed value of 
the oven-dried shrinkage a\ (= 0.2771%). In the present 
simulation, we assumed £ 3 m = £ 2 m = £ 3 m , and £ f x = e 2 for 
convenience. 

Optimized values of the increments in £, m , e[, and e 2 

were obtained by the simulation as displayed in Table 6. 
In our previous report, we succeeded in simulating the 
observed relationships between the longitudinal and the 
tangential swellings, and the MFA in the clear wood speci¬ 
men of sugi (Cryptomeria japonica ) by supposing that £ 3 m = 
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s 2 = 12%—15%, and e[ = e 2 = 0%-l%. 12 In the present 
simulation, optimized e[ and e 2 became very small but nega¬ 
tive, which means that the polysaccharide framework 
bundles in the SI and S2 layers tend to contract in the 
direction parallel to the cellulose molecular chains in spite 
of increasing moisture content in the cell wall. It is impos¬ 
sible for the authors to provide any comment on this strange 
result at this stage. However, their absolute values were 
very small compared with the increments in e™ and e™. 

We also simulated the oven-dried shrinkage of the 
G-fiber (a g L = 1.3026%), and optimized the value of 
the increment in e 3 . In this simulation, it is assumed that the 
increment in each of e™, e 2 , e[, and s 2 takes an identical 
value between the N-fiber and the G-fiber (see Subsidiary 
condition 3). For convenience, we assume £ : m = e 2 = £ 3 m in 
this simulation. Thereafter, we optimized the value of the 
increment in e 3 to obtain the observed value of a g L . Results 
are displayed in Table 7. The optimized value of the incre¬ 
ment in e 3 became a large positive value which is quite 
different from those in e\ and s 2 . 


Discussion 

Young’s modulus of the green G-layer (EJ) 

According to the results from Eqs. 7, 9, and 10, the pre¬ 
dicted Young’s modulus of the green G-layer (E<f) became 
2.15 times as large as that of the lignified layer in N-fiber 
(£n), and 1-74 times as large as the one in the G-fiber (EJ). 
In any case, we can say that the longitudinal Young’s modu¬ 
lus of the G-layer becomes larger than that of the lignified 
layer in the G-fiber and N-fiber. The predicted value of the 
Young’s modulus of the lignified layer in the G-fiber (E<J) 
became slightly larger than the one in the N-fiber (E^). This 
is because we calculated the value of in due consider¬ 
ation of an experimental fact that the MFA of the S2 layer 
in the G-fiber was a little smaller than in the N-fiber. This 
may be one of the factors to increase the Young’s modulus 
of the TW xylem. 

It is well known that the TW becomes very stiffer in the 
longitudinal direction as it dries. The increase of the 
Young’s modulus of the TW xylem due to drying is highly 
correlated with the percentage of the G-fiber in the fiber 
domain. 4 This suggests that the G-layer becomes abruptly 
rigid as the water molecule is released. However, the valid¬ 
ity of this suggestion remains to be confirmed. 

Strangely, the predicted value of the relative crystallinity 
in the framework bundle of the oriented polysaccharide in 
the G-layer was somewhat smaller than that in the second¬ 
ary wall (see results from Eq. 9). According to the formula 
derived in our previous report, 12 the Young’s modulus of 
the G-layer is highly dependent on the ratio of the cellulosic 
component. In the present simulation, we assumed it to be 
90% in the G-layer, which may be a little larger than the 
true value in the green G-layer. It is likely that the G-layer 
contains a substantial amount of noncrystalline polyose, 
e.g., hemicellulose. As another possibility, we note the fact 


that the green G-layer is highly swollen by water, which 
causes an apparent decrease in the relative crystallinity of 
the cellulose in the green G-layer. Hitherto, we have re¬ 
ferred to the classical data of Norberg and Meier 7 on the 
chemical and physical properties of the G-layer in aspen. 
However, we need to critically verify their conclusion for 
various species. 

Growth strain in the G-layer (e 3 ) 

The simulated value of e 3 is more negative than that in the 
lignified layer (s 1 and e f 2 ). This indicates that a large con¬ 
tractive internal strain originates in the polysaccharide 
framework of the G-layer in the direction of the cellulose 
molecular chain, which causes a high longitudinal tensile 
growth stress in the TW xylem. 

Shrinkage and swelling of the G-layer due to 
moisture adsorption 

Many authors have considered that the polysaccharide 
framework does not swell or shrink by the action of mois¬ 
ture adsorption. However, the present simulation shows 

r 

that the value of e 3 , which is the swelling ability of the 
polysaccharide framework in the G-layer, becomes a large 
positive value. Conversely, the polysaccharide framework 
in the G-layer tends to shrink in the direction parallel to 
the cellulose molecular chain during moisture adsorption. 
This means that the high longitudinal drying shrinkage in 
the TW xylem is induced by the drying shrinkage of the G- 
layer in its axial direction. Recently, Clair and Thibaut 28 
observed by scanning electron microscopy that the dried G- 
layer tends to be depressed from the surrounding lignified 
layer, which supports the predicted results in the present 
simulation. 
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Appendix A Deriving Eqs. 2 

We denote the number of the G-fiber in the wood 
fiber domain with an area A as G and that of the N-fiber 
as N, provided that G + N = X. We set the following 
assumption. 

Assumption A: The thickness of the lignified wall in the G- 
fiber is identical to that of the N-fiber regardless of the 
measuring position. 

This assumption is not so inappropriate to the wood fiber 
domain in the real xylem because the observed values of s 


and XIA (= /) become almost unchanged regardless of 
measuring positions, as seen in Table 1. Moreover, we set 
the following assumptions. 

Assumption B: The diameter of the G-fiber is similar as that 
of the N-fiber. 

Assumption C: Cellular arrangement in the crosscut surface 
of the xylem takes a tessellation structure consisting of a 
polygonal cell. 

Then, we can connect the area ratio of the lignified layer [s], 
and that of the gelatinous layer [g] to p 0 , p 1? p 2 , and p 3 under 
assumptions A, B, and C. 

It may be a little hasty to apply calculated values of p 0 , p 1? 
p 2 , and p 3 to the simulation using Eqs. 1 because the crosscut 
shape of the G-fiber model displayed in Fig. 1 is close to 
circular. However, we know that the hexagon is the most 
closely allied to the circle in shape among the polygons that 
constitute the tessellation arrangement. Then, we set the 
following assumption. 

Assumption D: Crosscut shape of the wood fiber in the 
cellular arrangement is a hexagon with an area of 

(3V3/2)r 0 2 as displayed in Fig. 2. 

We denote the thickness of the lignified layer as 
|V3/2j(r 0 — r 3 ), and that of the G-layer as (S/2]{r 3 - r 4 ). 
Distances from the central point of the hexagon to 
the lignified and gelatinous layers are denoted as |V3/2jr 3 

and |V3/2jr 4 , respectively, s can be given as in the following 
formula: 



Fig. 2. A hexagon-shaped model of the wood fiber used for calculating 
Eqs. 2. T, tangential direction; R, radial direction 





























5 = |V3(r 0 2 - rl)x/A = 1 - |>/3 • r 3 2 • /. (Al) 

In a similar manner, g is given as the following formula: 

g = |V3(r 3 2 - r 4 2 )G/A = (l - - |V3 • r 4 2 • j. (A2) 

XIA, GIA , s, and g can be determined experimentally as 
displayed in Table 1. From Eqs. Al and A2, we can obtain 
r 3 and r A as follows: 

r > = | 1 j 2 '/3( 1 “ *)y. 

r ‘ ■ I A 1 ~ s) i ~ 4 


■ (G # 0). 


(A3) 


If we denote r 0 /r 1? rjr 2 , r 2 /r 3 , and r 3 lr A as p 0 , p 1? p 2 , and p 3 , 
respectively, we obtain the following equation: 


P 0 P 1 P 2 


VT^ 


p 3 


( 1 -,) 


W ~s)-g-f/N, 


(for A g # 0), p 3 = 1 (for A g = 0), 


(A4) 


where / = A/A, and A = G/A. 


Appendix B Integration (Eq. l'J is not affected by the 
functional shapes of t-dependent variables S lf S 2 , S 3 , 

8™, 82, £ m , 8^ 8^, and 83 

We introduce functions cp A , cp 2 , and cp 3 which vary from 0 to 
1 in the range of 0 < t < T 3 as follows: 



(0 < f < T x ) 

(T\ — t — t 3 )’ 



(0 < t < r 2 ) 
(r 2 < r < r 3 )’ 


0 (0<r<7;) 

<2(r) (7) < r < T 2 ), 

1 (r 2 <r< T 3 ) 

(Bl) 


where E(0, 2(0, and P(0 are monotonously increasing and 
differentiable functions, which vary smoothly from 0 to 1 in 
the respective domains. With reference to Subsidiary condi¬ 
tions 1 and 2, we assume the following condition as the 
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functional shapes of Fdependent variables S t (t), £ z m (0, and 
0(0 (/ = 1, 2, 3): 


Si{t) = ki ■ tp^t), £,“(?) = m, ■ <Pi(t), e, f (f) = n, ■ (Pi(t). 
(i = 1, 2, 3) 


(B2) 


where k h m h and n { are constants. It is reasonable to assume 
this condition if these Fdependent variables change the val¬ 
ues smoothly during maturation of the matrix skeleton in 
their respective integration intervals. Then, by substituting 
Eq. B2 into Eq. 1', we obtain the following expression: 


dt 




+ 


+ 


3 

J 

J 


dt 


dt. 


(B3) 


where 




gi 


WO) = 


£3 wo) 


+ «l/l4(p)| 5l = tl9l(f ) (0^ 

^2 fl 2 (p)| S 2 = ^2 • 

+ « 2 /l 5 (p) 
m 3 ' /u(p) 

+ ^/i«(p)| s ,-* 3 ..„(,) (r 2 <f<r 3 ) 


S 2 = k . 2 -(fiil'l dl 

S 3 =k 3 -q> 3 (t) 


T t ) 


r 2 ) 


From Eq. Bl, Eq. B3 is modified into the following 
expression. 


A = G g\P(t)}dP(t) + l T ; g2 [Q(t)\dQ(t) 

T (B4) 

A->[*( 0 M 0 - 

P(t), 2(0. and i?(0 change the values from 0 to 1 monoto¬ 
nously and continuously for elapsed time t in their respec¬ 
tive integration intervals. We can rewrite Eq. B4 as the 
following expression: 



This result indicates that the integration value in Eq. B5 
does not depend on the concrete values of 7\, T 2 and T 3 . 
Furthermore, it is not affected by the functional shapes of 
the Fdependent variables S l9 S 2 , S 3 , 8™, £™, e[, £ f 2 , and s f 3 if 

we assume the condition described by Eq. B2. 





















